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Abstract 

In this paper we calculate the Rayleigh-Brillouin spectrum for a relativistic simple fluid 
according to three different versions available for a relativistic approach to non-equilibrium 
thermodynamics. An outcome of these calculations is that Eckart's version predicts that 
such spectrum does not exist. This provides an argument to question its validity. The 
remaining two results, which differ one from another, do provide a finite form for such spec- 
trum. This raises the rather intriguing question as to which of the two theories is a better 
candidate to be taken as a possible version of relativistic non-equilibrium thermodynamics. 
The answer will clearly require deeper examination of this problem. 
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I. INTRODUCTION 



It is a well known fact that light scattering by a simple fluid in equilibrium at 
a certain temperature T and pressure p is one conclusive test to verify Onsager's 
linear regression of fluctuations hypothesis a basic assumption in classical irre- 
versible thermodynamics j3, S|. Not only it constitutes the core behind the proof of 
Onsager's reciprocity theorem but it also guarantees that the equilibrium state of 
the fluid is stable under such fluctuations. Thus, it is legitimate to ask if the various 
formulations of irreversible relativistic thermodynamics so far available are at grips 
with such hypothesis. Even in an indirect way, this would indicate that such theories 
can be tested experimentally. Indeed, in the classical case one measures the so-called 
dynamic structure factor of the fluid S (g, u) which represents the energy scattered 
by the fluid from an incoming wave of wavelength A, as a function of frequency. The 
outcome of this measurement is the well known Rayleigh-Brillouin (RB) spectrum 

nriR 

[4, 15|, 16|. Its central peak, the Rayleigh peak, has a width proportional to the ther- 
mal diffusivity Dt = n/poCv where k is the thermal conductivity, po the equilibrium 
density and Cy the specific heat at constant volume. This peak represents the in- 
tensity of the thermal (entropy) fluctuations. Symmetrically located with respect to 
this peak there appear two peaks, the Brillouin peaks, which represent the fluctu- 
ations arising from the mechanical dissipative processes of the fluid, the sound or 
light absorption. They are located at u = ±Cok from the central peak, Co being the 
velocity of sound, if the probe is a sound wave. Their width is given by the famous 
Stokes-Kirchhoff's formula namely. 



where 7 = Cp/Cy, f] and ^ being the shear and bulk viscosities, respectively. 

The important feature here is that the precise form of this spectrum can be 




(1) 
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obtained by solving the linearized Navier-Stokes- Fourier equations of hydrodynamics 
for the perturbations (or fluctuations) (5T, 5p and 5u present in the fluid due to its 
microscopic structure. 

Now, it may be so that this experiment per se could not be easily carried out in a 
laboratory in a relativistic regime for technological reasons. However, one surely must 
expect that the corresponding linearized equations of relativistic hydrodynamics lead 
to a relativistically modified spectrum which reduces to its classical counterpart in 
the non-relativistic limit. This is precisely the motivation of this paper. We wish to 
calculate the RB spectrum for the linearized relativistic hydrodynamic equations that 
arise in three cases: using the Eckart-Landau Lifshitz formalism p, 18|, a relativistic 
generalization of Meixner's theory jo] and considering the equations obtained by the 



authors 
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ll| when the acceleration term in Eckart's theory is expressed in terms 



of Vp using Euler's equations. We will refer to this case as the modified Eckart's 
theory. 

In Sects. II-V we shall establish the system of linearized relativistic fiuid equations 
for the three alternatives mentioned above and analyze the modifications to the RB 
spectrum in each case. Section VI is devoted to the discussion of the results and 
final remarks. 



II. MEIXNER-LIKE FORMALISM 

The first formalism we wish to analyze consists of the relativistic generalization 
of Meixner's formalism Q|. In it, the heat flux is not considered as part of the 
momentum-energy tensor but is included in a separate total energy flux conserva- 
tion equation. As a result, the constitutive equation for the heat flux retains its 
Fourier-type structure. The hydrodynamic equations for this formalism have been 
obtained elsewhere j^, The fluctuations, here denoted by a 5 preflx, evolve to 
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the equilibrium state following the linearized version of such equations. Thus, the 
dynamics of the fluctuations is given by 



Sn + no66 = 
Po59 + -^VHn + ^V^5T - AV^59 = 

6f + 56 - DtV^5T = 



(2) 
(3) 

(4) 

The second equation is a balance equation for the longitudinal component of the 
fluctuations in the hydrodynamic velocity given by 66 = du':'^,, where u'^ is the hydro- 
dynamic velocity four- vector. This equation is obtained by calculating the divergence 
of the momentum balance equation, a procedure which decouples the transverse mode 
whose dynamics has been already analyzed in a separate work jl^. Here n is the 
particle number density, T the temperature, kt the isothermal compressibility, (3 the 
thermal expansion coefficient, C„ the heat capacity at constant particle density and 
A = ( + Arj/S. We have defined po = (^o^o + Po) 1^ where and po are the internal 
energy and pressure respectively. The naught subscripts denote equilibrium values, 
the semi-colon a covariant derivative and a colon a component of a gradient. Greek 
indices run from 1 to 4 and Latin ones from 1 to 3. 

As mentioned above, the procedure to obtain the spectrum is the standard one 
and involves calculating the dispersion relation arising from the determinant of the 
Fourier-Laplace transformed hydrodynamic system of equations. For Eqs. ©-([l]) 
we obtain that 



no 



■^^q^ s + 
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s + Drq^ 



(5) 
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which yields a cubic dispersion relation which can be written as follows 

+ a2S^q^ + s [a^q^ + a^q^^ + asg^ = (6) 
where the coefficients are given by 

A ^ 
a2 = — + Dt 
Po 

as = ^Dt (7) 
Po 

7 
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as 



ktPo 
Dt 



ktPo 

and we have used the relation ^ = = 7 — 1. One can easily show that Eq. 
([3]) has one real root given by 



and a pair of conjugate roots 

■52,3 



'^5 2 io\ 

Si = q 8 

a4 



The analysis follows exactly as in the non-relativistic case j^-Q] where it is shown 
that Eqs. ([8]) and ([9]) are valid up to terms of order g"^. Recalling that S (g, lS) is the 
density-density self correlation function, we may plot the ratio between the dynamic 
and static structure factors, Si^q^iS) / S (q) as a function of u for a fixed q which 
should yield three peaks given by the roots of the cubic equation. The mean width 
of the central peak, the Rayleigh peak, is determined by the real root given in Eq. 
([Hj) and thus, for the Meixner, case we obtain a width 

^RM = — (10) 

7 
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A correction due to the modified value of the thermal conductivity k for relativistic 
fluids will arise and thus one expects to observe a change in the width of the peak. 

The location and width of the other two peaks, the Brillouin peaks, are determined 
by the conjugate roots given in Eq. Q. The location of the peaks is given by the 
imaginary part while the width is given by the real part. In this case the doublet 
appears at = i^/olg so that from Eqs. ((71) it follows that 



= ±J—^q (11) 

V i^tPo 

One should find a shift in this location due to the relativistic value of po- The width 
of the doublet is given by the real part of S2,3, that is 

(^2 0.5 \ „,2 



and in this case we have 



Asm = l(^ + ^^Dt] (13) 

where once again, a correction due to the relativistic values of Dt and po is expected. 
In both Eqs. ffTOll and (fT3ll one recovers the non-relativistic expressions when c — > cxd 
as in this case Po ^ Po, the equilibrium density. 



III. ECKART'S FRAMEWORK 

Eckart's theory for relativistic fluids [3| is based on the construction of an 
momentum-energy tensor where heat flux is included. As a consequence, in or- 
der to satisfy the second law of thermodynamics, he proposed a rather controversial 
constitutive equation for the heat flux in which a hydrodynamic acceleration term 
is included. This proposal has been claimed to render the theory unphysical and 
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motivated the use of extended theories as alternatives [l3|-[l5|, as has been thor- 
oughly discussed 1^, [ill, [l^. The aim of this section of to present the effect of such 
constitutive equation in the structure of the RB spectrum. Thus, once again, the 
starting point is the linearized set of equations for the fluctuations in a relativistic 

n 

fluid, now within Eckart's theory, which are shown in Ref. pj| to read as 



6h + no59 = 



(14) 



6f + -^^69 

noCnKT 



2.^ DtTo^^ 







(15) 



(16) 



Indeed, the two terms -^69 and -^^69 in Eqs. ^ and ^ come from Eckart's 
proposal for the heat flux constitutive equation depending on the hydrodynamic 
acceleration through the term —■^u'^ (see Eq. (fTOll ). All quantities appearing in 
these equations are the same ones that appear in Eqs. (EMD- 

Proceeding as in the previous case, we analyze the dispersion relation which is 
now given by 



s no 

- s + DtQ^ 

which yields a quartic polynomial, namely 

bis^ + s^ + b2S^q^ + s (bsq"^ + a^q"^) + b^q^ = 

with the coefficients given by 



(17) 



(18) 
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, ADt 



b. 



'4 



Po 

7 

h = — - 

ktPo 

We now attempt a rather intuitive but accurate solution to Eq. (fTSll . Notice that 
the coefficient of s"^ is very small. Neglecting it, we can readily identify three roots, 
namely, 

^5 2 

Si = -—q 
64 

Now, we assume that these three roots are still approximate solutions to the quartic 
and find the fourth root by using the property that, since the coefficient of is equal 
to one, Z]f=i Si = — ^, and thus 



c Po , 
S4 ^ ^ + 



Po \ C^Pol^T, 



Since [3 / kt < 0, the fourth root S4 is always positive. A real positive root in the 
dispersion relation yields an exponential growth in the structure factor instead of 
a finite spectrum. This behavior is unphysical and simply implies that the RB 
spectrum does not exist in Eckart's formalism, even in the non-relativistic limit. 

IV. MODIFIED ECKART'S THEORY 

The system of equations we consider in this section is obtained, as in the previous 
one, from a momentum-energy tensor which includes relativistic heat fiux terms. 
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The key difference fiere is tliat tlie constitutive equation introduced for tlie heat flux 
is obtained through the following argument. According to Eckart, such constitutive 
equation is given by 

Jm = -^k{j^' + 1^'^) (19) 

Iq violates the tenets of classical irreversible 



where the second term, as argued before 



thermodynamics since it is neither a thermodynamic force nor a flux. Further, it has 
another serious drawback namely, it raises ii'^ to the category of a state variable, a set 
already chosen to be given by n, u'^ and T. Thus, the resulting set of hydrodynamic 
equations would be overdetermined. Therefore, to keep Eq. ffTOl) to first order in the 
gradients, we eliminate it'^ using Euler's equation 

Pou" = -p,,h^'' (20) 

Now, according to the local equilibrium assumption, 

P,^, = —T^ + -^n^^ (21) 

and thus, substitution of Eq. (i2Ti ) in Eq. ( flQl l yields 

4l = - fl + ^l) - (22) 



or 



Jfg] = -LttT'' - LnTw' (23) 
where Ltt is an "effective thermal conductivity" given by 

Ltt = n(l + -^^] (24) 

and LnT a new transport coefficient given by — ^"^"o- which has no classical counter- 
part. We would like to remark that an equation similar in structure to Eq. f l23l l was 



already derived by Landau and Lifshitz 
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a. 



Clearly these transport coefficients can also be calculated from a kinetic model 
but we shall discuss such calculation in a separate paper. The equations in this 
formalism need not be justified in detail since they follow simply from the method 
used to obtain Eckart's equations incorporating the terms arising from Eq. f l23ll for 
the heat flux. Thus we obtain that 



Sn + no5e = 



(25) 



Po6e + -^V^Sn + ^W^ST 



L 



TT 



6f + -^59 
and give rise to the dispersion relation, 



^V'5n = 
noCn 



(26) 
(27) 



no 



noCn ^ 



To/3 

noCnKT 



which can be written as 



(28) 



+ d2S^q^ + s (d^q'^ + d^q^) + d^q^ 



(29) 



where 



A Ltt a [^Tq \ TlQLnT 

02 = — H I 1 



di 



d^ 



1 



ALfj- 



KTPonoC„ 



riQCnPo 

- ^ 
ktPo 

{Ltt - PnoLnT) 



(30) 

(31) 
(32) 

(33) 
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Once more, following the steps outlined for the two previous cases, one can identify 
the modifications to the spectrum, as before, by analyzing the roots. For this case 
the width of the Rayleigh peak is given by 

^RS = ^q' = {Ltt - /3noLnT) (34) 

The shift in the Brillouin doublet is the same as in the Meixner case. In the particular 



case of an ideal gas, the properties of po 



17| guarantee that the position of the peaks 



will never exceed ±cg. On the other hand, the width is significantly modified. We 
now obtain 

2 [po noCn\ 7 c^i^tPq) \Cnl c^PoJ) 

Equations (l34ll and fl35l) deserve further attention. The former predicts a modi- 
fication in Rayleigh's peak which changes its width due to the effective thermal 
conductivity given by Eq. (l24ll and the presence of L^^; both, as stressed above, are 
strictly relativistic effects. This is quite different from Meixner's like theory where 
the correction arises only from the relativistic value of k/c„. Moreover, the shape 
of Brillouin's peaks is further altered due to several relativistic terms as can be seen 
from Eq. (I35|). 



V. SUMMARY AND FINAL REMARKS 



In the previous section, the modifications to the RB spectrum according to the 
three versions of relativistic irreversible thermodynamics have been explored. The 
difference between the modified Eckart's theory result analyzed in Sect. IV and 
Meixner's case analyzed in Sect. II, should be emphasized. The latter one does 
not have a density gradient in "Fourier's equation" which, as shown in Eq. ( l22l) . is 
strictly a relativistic factor. This poses an intriguing question namely, in both cases 
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which are alternative versions of a relativistic non-equilibrium theory, we predict the 
existence of a RB spectrum. These spectra are different in both cases, but contain 
modifications in comparison to the classical spectrum and both reduce to it in non- 
relativistic limit. Which theory is the correct one? If we believe in fundamentals we 
would be inclined to choose the version that is consistent with the results obtained 
from kinetic theory and thus, as we have shown in previous work [l^, the second 
theory prevails. However, Eckart's modified theory still contains the heat flux as a 
component of the momentum-energy tensor which still is debatable. On the other 
hand if we wish a phenomenological theory that contains a density gradient in the 
heat flux within the Meixner's like formalism we would face the problem of how 
to introduce such term. The answers to these puzzles are still open. We feel that 
Eckart's original approach may be discarded but the flnal answer as to which is 
the appropriate version of a relativistic non-equilibrium thermodynamics is still a 
challenge both theoretically and experimentally. 

To facilitate the various results available for the RB spectrum we are summarizing 
them in the Appendix. 

APPENDIX 



The modifications to the RB spectrum in the three frameworks considered are 
summarized in the following table 





Rayleigh's peak width 


Brillouin peaks shift 


Brillouin peaks width 


Non-Relativistic 


1 ^ 


±^/ ^ 1 




q^ 


Meixner 


2t 2 
7 ^ 


±,/-T q 


^ {i + ^^-) 




Eckart 


No spectrum 


No spectrum 


No spectrum 


Modified Eckart 




±J^'' q 
Y PoKr ^ 


^ ( A Ltt ( 7-1 .^10^] L 




^l^^ {Ltt PnoL^r) 





where calligraphic fonts are being used for relativistic transport coefficients in order 
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to distinguish them from the non-relativistic ones. 
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